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METHO D OF OPTIMIZING PAR AMETER VALUES IN A PROCESS OF 
PRODUCING A~PRODTJCT 

«'? =— 

BACKGROUND OF THE INVENTION 
Field of the invention 

The invention relates to the process optimization field, and more 
particularly to a method of optimizing parameter values in a process of 
producing a product which is characterized by properties affected by the 
selected parameter values. This invention is applicable in different 
industries, such as the pharmaceutical, chemical, cosmetics, plastics, 
petrochemical, agriculture, metallurgy and food industries, as well as many 
other commercial and industrial applications. 
Description of prior art 

Processes for production of complex compositions such as those 
found in many pharmaceutical products generally require the mixing of 
many ingredients according to specific process parameters regarding 
formulation and production technology, to provide the product with 
properties at a level offering satisfactory performance according to 
predetermined specifications. In such complex production processes, it is 
not unusual that some process parameters involved exhibit interfering 
effects on the desired properties, further complicating the process design. 
Where possible, the designer may try to adapt the set of process 
parameters from known data derived from previous similar processes, 
and/or rely on conventional trial-and-error experimental schemes to 
optimize the set of process parameters values, in order to meet the product 
specifications. However, as the processes become more complex, 
optimization in such multidimensional space with high accuracy 
requirements turns out to be an extremely difficult task, even for the highly 
skilled designer. That limitation is particularly problematic in the design of 
pharmaceutical products, where one or more active substances mixed with 
a variety of excipients (e.g. carriers) must be produced in the form of a 
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stable and highly effective standard delivery system such as a tablet, 
capsule, suspension, cream or injection, or even controlled release 
systems such as skin carriers and implants. 

In the past years, many techniques have been developed to assist 
5 the process designer or formulator in optimizing values of parameters 
governing processes. These techniques aim at quantify existing relations 
between parameters and associated desired product performance 
characteristics. A conventional technique known as the Full Factorial Matrix 
(FFM) method consists of statistically deriving a behavior relations for the 

10 properties from a set of experimental runs of the process using selected 
initial values for the parameters. The established model being generally 
nonlinear, optimized parameter values are then derived using an 
optimization method such as the Multisimplex method described in 
"Practical Methods of Optimization" J, Wiley & Sons, Chichester, 2d, 

15 (1987), which essentially consists of linearizing the behavior functions 
related to the parameters according to straight lines or planes of different 
random directions. For any given property behavior relation of n 
parameters to be optimized in order to either minimize or maximize that 
behavior relation with or without constraints on the parameter values, a 

20 recursive estimation of the property is then performed using an initial set of 
parameter values according to a selected direction, until the obtained value 
for the property does not significantly vary in that direction. Then, a last 
unfavorable set of parameters is used as a new starting point for a 
following recursive estimation according to a different direction. Successive 

25 recursive estimation steps are performed until the resulting value for the 
property no longer significantly vary in any new direction. When applied to 
a model comprising a plurality of property behavior relations, the 
Multisimplex method allows a unique objective function to be created by 
proper transformation of the relations to adapt to different scales and/or 

30 units and by associating a relative importance weight to each property, 
either subjectively or through fuzzy logic algorithms. 
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The known optimization processes based on Full Factorial Matrix- 
Multisimplex methods suffer from several drawbacks. As a general rule, the 
number of experimental runs required to obtain a model of sufficient 
reliability is proportional to the total number of significant parameters 
5 involved. Therefore, the cost and time frame of the experimental work will 
therefore be essentially proportional to the number of runs required. 
Although a variant of the method known as the Fractional Factorial Matrix 
has been proposed in order to reduce the number of runs to be performed, 
the provided reduction of experimental runs may not significantly reduce 

10 the total cost and time frame of the work required to complete the design of 
a complex product involving many production technologies. While adequate 
formulations complying with constraints imposed on the parameter values 
can nevertheless be obtained, these formulations generally cannot be 
qualified as optimal when comparing actual property performance with 

1 5 desired property values set forth in the product specifications. 

A technique which attempts to improve parameter optimization in 
process design is disclosed in European Patent Office laid-open patent 
application publication number 0,430,753 dated June 5, 1991 and in US 
patent No. 5,218,526 issued on June 8, 1993 to Mozzo. According to the 

20 technique in Mozzo, from a set of property relations expressed in terms of 
parameters which is obtained by standard statistical methods using the 
results of a number of experimental runs of the process, a corresponding 
set of property relations expressed in terms of weighted parameters is 
derived. For each actual value of a parameter, a first weighting is 

25 expressed as the ratio of: (a) the deviation of the actual value from the 
mean value of the parameter over the experimental range, on (b) the range 
between extreme values for that parameter over the experimental range. 
Then, a goal function is established in term of deviations between weighted 
values of property values as estimated by the property relations and 

30 corresponding weighted values of specified goal values for the properties. 
For each goal value of a property, a second weighting is expressed as the 
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ratio of: (a) the deviation of the actual value from the mean value of the 
property over the experimental range, on (b) the range between extreme 
values for that property over the experimental range. Then, according to a 
recursive geometric algorithm aimed at successively minimizing the 
5 established goal function, a set of optimal parameter values is generated. 
While being an improvement over the conventional Full/Fractional Factorial 
Matrix - Multisimplex methods regarding the capability to consider specified 
goal values for the properties, the weightings as taught by Mozzo do not 
reflect the relative importance of the properties involved, and that limitation 
10 may therefore affect the convergence of the algorithm toward an optimal 
solution. 

A review of modern techniques and software systems for the design 
of pharmaceutical product formulations is given in "Intelligent Software 
System For Pharmaceutical Product Formulation" R.C. Rowe, 

15 Pharmaceutical Technology, March 1997. In that paper, expert systems, 
rule induction algorithms, case-based reasoning algorithms, neural 
networks and genetic networks are presented as modern tools for 
supporting formulation design, and a number of available software systems 
using some of these tools are summarized. As indicated in the Rowe 

20 paper, although a knowledgeable expert system could be a powerful tool to 
assist the process designer in the formulation task, its development is 
generally a high risk, time consuming and expensive process. Rule 
induction is a knowledge-based algorithm which allows hierarchical 
classification of objects, using statistical methods which are found generally 

25 effective only if the input data is continuous, which is often not the case in 
practice. Moreover, since rule induction is limited to establishing whether or 
not a given object is close to another, it generally cannot provide an optimal 
solution. Case-based reasoning is a knowledge-based iterative technique 
which can be used to design formulations, which consists of matching the 

30 desired specifications for the product with the specifications of the most 
relevant known formulation(s), and adapting the selected formulation(s) as 
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necessary, followed by an evaluation. Although effective for optimizing the 
parameters of a variant process from a family of similar processes and 
corresponding formulations, case-based reasoning generally cannot be 
used where the design of a significantly different formulation is 
5 contemplated. As to neural networks, in which each neuron input is 
modified by a weight associated with that neuron, they appear to be 
effective tools for assisting formulation design only in cases where no 
constraint applies on either the parameter or property values, such cases 
being rarely found in practice. Finally, regarding the genetics algorithms, 
10 they are cyclic methods based on Markov chains for predicting from a 
starting point a solution likely to result from a sequence of operations, in 
order to allow making changes to obtain a desired solution. Since these 
changes are generally made arbitrarily, in most cases, the resulting solution 
cannot be considered as optimal. 

15 

SUMMARY OF THE INVENTION 

It is therefore an object of the present invention to provide a 
systematic method of optimizing parameter values in a process for 
producing a product which minimizes the number of experimental runs 
20 required to obtain an optimal solution complying with the product 
specifications. 

According to the above object, from a broad aspect of the present 
invention, there is provided a method of optimizing parameter values in a 
process of producing a product, the process being essentially controlled by 
25 a set of n parameters X i affecting a set of k properties Y } 

characterizing the product. The method comprises the steps of: i) 
assigning values to a set of k property weights Wj representing relative 

importance of the properties Y } for the characterization of the product; ii) 
establishing property behavior mathematical relations giving an estimated 
30 property Ye } for each property Y } in terms of the parameters X { from given 
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parameter data and associated property data; iii) using the property 
weights w ; to establish a goal function in terms of property weighted 

deviations between the estimated properties Ye } and corresponding 

specified goal values for the properties Y } \ and iv) optimizing the goal 

5 function to generate a set of n optimal parameter values for the 
parameters X r 

According to a further broad aspect of the present invention, there is 
provided a method of producing a pharmaceutical product using optimized 
process parameter values, the process being essentially controlled by a set 
10 of n parameters X i characterizing a formulation for the product, the 
parameters X i affecting a set of k properties Y } characterizing the 

product. The method comprises the steps of: a) conducting a number of / 
of experimental runs of the process each using a selected distinct set of 
values for the parameters X i covering substantially all extreme values 

15 within a chosen range of values for each one of the parameters X n 
wherein / is at least equal to n + 1 and is substantially less than a number 
used in the Fractional Factorial Matrix method; b) measuring values for the 
properties Y } characterizing the product in each of the / experimental runs, 
whereby parameter data and associated property data are obtained from 

20 the selected distinct set of values for the parameters X { and the measured 
values for the properties Y } , respectively; c) determining an importance of 
the properties Y } for the characterization of the product, comparing the 
importance of the properties Y. relative to one another, and assigning 
values to a set of k property weights Wj representing a relative importance 

25 of the properties Y } for the characterization of the product; d)calculating a 
set of optimal parameter values for the parameters X. t using the measured 
values for the properties Y. and the assigned values of the set of k 
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property weights w y . ; and e) producing the pharmaceutical product using 
the optimized process parameter values X. calculated in the previous step. 

BRIEF DESCRIPTION OF THE DRAWING 

5 The invention will be better understood by way of the following 

detailed description of a preferred embodiment with reference to the 
appended drawings, in which: 

Fig. 1 is a block diagram of a software system that can be used to 
carry out the method according to the present invention according to the 
10 preferred embodiment; and 

Fig. 2 is a flow chart representing the preferred embodiment of the 
method according to the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

15 In the following description, a preferred embodiment of the present 

invention applied to product formulation design will be described. However, 
it is to be understood that the present invention can be also be used to 
optimize parameter values of processes related to the production of many 
types of products which cannot be associated with a formulation, while 

20 being characterized by a number of properties affected by process 
parameters, such as biotechnological products, electronic components, etc. 

Referring now to Fig. 1, there is generally designated at 10 a 
computer system which is programmed to carry out a method according to 
the present invention. The system preferably comprises a knowledge base 

25 12 where prior formulation/process data and competing products data are 
stored. For the purpose of pharmaceutical formulation design, knowledge 
base 12 contains process data related to ingredient proportions, 
experimental conditions and results over time, production technologies 
used, etc. The system 10 further comprises a property weighting module 14 

30 which generates a weight value for each one of a number k of identified 
properties according to an initial modeling of the problem and property 
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comparison data presented to the module 14. System 10 further comprises 
an evaluation module 18 fed by the property weights generated by module 
14, to generate a global relative importance vector of dimension [k] for the 
k properties. System 10 is provided with an experimental data entry 
5 module 16 through which property values obtained from experimental runs 
using different sets of parameter values for the process can be entered and 
stored for later use by several modules of system 10. Linked to receive 
data from modules 12, 14 and 16, is an evaluation module 18 which can 
generates a ranking of the sets of property values selected from the 

10 knowledge base and the optimal set of parameter values obtained through 
optimization. System 10 further comprises a parameters reduction module 
22 to retain only those parameters having a significant effect on the 
considered properties. Module 22 is particularly useful in cases where the 
number of parameters involved is relatively large, usually greater than 8 

15 where a computer provided with a standard high-performance 
microprocessor is used. The S-Plus™ statistical software from MathSoft 
may be used in module 22 to carry out the Stepwise method to select the 
variables . System 10 is further provided with a parameter interaction 
module 20, the function of which consists of identifying by statistical 

20 analysis from experimental data, which remaining parameters are 
significantly correlated. The S-Plus statistical software from MathSoft can 
also be used to program module 20 in which the appropriated correlation 
methods are applied to the data. It is to be understood that module 20 is 
unnecessary where all parameters are independent one another. 

25 Modules 16, 20 and 22 are linked to a property behavior models 

mocjule 24 that uses experimental data, parameter interaction data and 
remaining^ignificant parameters for determining an optimal mathematical 
model for each properf^^i^hjs likely to better estimate that property. The 
model data as generated at module*244s^fed to a property behavior relation 

30 module 26 that also receives experimentaT^ta^rom module 16 to 
statistically estimate polynomial coefficients to be incorpoi^ted-Ayithin the 
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established property behavior models, thereby generating a behavior 
relation for each property. The S-Plus statistical software from MathSoft 
rnay^bei used to program module 26 to apply the appropriate regression 
methods tolfife^ata^System 10 is further provided with a goal function 
5 module 28 linked to property^Weighting module 14 and property behavior 
relation module 26 to generate, from specified goal values for the 
properties, a goal function in term of property we^ghteB^tieviations between 
properties as estimated by the behavior relations and the ^orre^onding 
specified goal values for these properties. 

10 An optimization module 30 is provided to optimize the goal function 

as established by module 28 through successive iterations and according 
to the type of each variable (discrete or continuous) and according to one 
or more ranges specified as constraints imposed on one or more optimal 
parameter values. Module 30 can be programmed using Matlab™ software 

15 supplied by The Math Works Inc to implement network optimization 
methods. Optimization module 30 is linked to the experimental data entry 
module 16 to transfer thereto the generated set of optimal parameter 
values, which module 16 also stores the actual property values obtained 
from an experimental run based on the set of optimal parameter values. All 

20 experimental data is then transferred to the evaluation module 18 as 
mentioned before. 

A preferred embodiment of an optimization method according to the 
present invention will now be described with reference to Figs. 1 and 2. As 
illustrated in the general flow chart shown in Fig. 2, the method comprises a 

25 first step 40 of assigning values to a set of k property weights w y 

representing relative importance of the k properties Y } for the 

characterization of the product, which properties are likely to be affected by 
the parameters of the process, from a modeling of the problem expressed 
as a hierarchical tree of these properties. Initial modeling and weight value 
30 generation are preferably performed using a method known as analytic 
hierarchy process (AHP), which was first proposed by T.W. Saaty, and 
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more recently described in "Using The Analytic Hierarchy Process For 
Decision Making In Engineering Applications: Some Challenges" 
Triantaphyllou et al, International Journal of Industrial Engineering: 
Application and Practice, Vol. 2, No. 1, pp. 35-44, (1995), which is 
incorporated herein by reference. 

The AHP method consists of building a hierarchical tree from all 
prSjpecties, with one or more hierarchical levels depending on existing 
relations betwfceqtoe properties. For each level, a pair-wise comparison 
matrix is built between th^foparties of this level and presented at an input 
10 of the parameters weighting module I^TsfiDwqJn Fig. 2, which executes in 
step 40. For each pair-wise comparison, the normalized eigenvector is 
derived associated with the higher eigenvalue. The compbn^nts of this 
eigenvector give the relative importance of each property called tffe^lQcal 
weight. Finally, the above normalized vectors are combined to find the 
1 5 global weight for each property. 

In a parallel direction, each pair-wise comparison is associated with 
a consistency index reflecting the transitivity relation between all 
comparison by pairs given by the formulator. Multi-criteria analysis software 
which is commercially available, such as Expertchoice™, Criterium™ or 
20 Ergo™, may be used to program module 14. For example, to one or more 
m main properties classified at a first (higher) level, may correspond one or 
more groups of properties classified at a second (lower) level, the latter 
properties being therefore identified as sub-properties. For each main 
property associated with a group of p sub-properties, a matrix of 
25 dimension [/? + lx/? + i] is built and filled, as a result of a pair-wise 
comparison between each property and sub-property, using relative 
importance values selected from a standard AHP scale. Next, a suitable 
algorithm performed by parameter weighting module 14 consists of first 
calculating the higher eigenvalue of the resulting numerical matrix, and 
30 then deriving a normalized relative importance vector of dimension [p + 1] 
by an estimation of the left principal eigenvector of that matrix associated 
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with the calculated main eigenvalue of the input matrix. The above 
algorithm is then applied to compare the m main properties of the higher 
level, from a pair-wise comparison matrix of dimension [/wxm] from which 
a normalized relative importance vector of dimension [m] is derived. 
5 Finally, the above normalized vectors are combined according to the 
hierarchical relations to generate a global relative importance weight vector 
for the k properties of dimension [m + or [k]. In practice, it is 

generally appropriate to retain only each group of sub-properties without 
the corresponding main property, the sum of the weights related to the 

10 retained k properties/sub-properties being always equal to unity. 

According to the next step, namely step 42, parameter data and 
property data values are provided, which data is obtained from 
experimental runs using different sets of parameter values for the process, 
the various values for each parameter being preferably selected according 

15 to an expected operation range within which an optimal parameter value is 
likely to be found. The parameters X t used in the experimental runs should 

cover the extremes of the expected operational range for each parameter. 
Generally, the number of formulation combinations required to determine 
an optimal formulation depends on many factors among which the more 

20 important ones are: 1) the formulation designer experience; 2) complexity 
of the formulation; 3) the availability of literature and experimental data 
available on the desired product; and 4) the analytical laboratory workload 
and throughput. According to the method of the present invention, the 
minimal number of experimental runs / to perform has been found to be 

25 equal to n + 1 , wherein n is the number of relevant parameters involved. A 
greater number of runs is certainly possible. Step 42 is performed by 
experimental data entry module 16 shown in Fig. 1. 

The method then comprises a step 44 of establishing property 
behavior mathematical relations linking the properties with the parameters 

30 and interactions thereof, in polynomial form. These property behavior 
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relations provide an estimated property Ye } for each of the k properties Y l 
in terms of a number n of parameters X i from the parameter data and 

associated property data provided at step 42. Step 44 is typically 
comprised of four sub-steps, namely 1) a parameters reduction step 
5 performed by module 22, 2) a parameters interaction analysis step 
performed by module 20, 3) a property behavior modeling step performed 
by module 24, and 4) a property behavior relations generating step 
performed by module 26, as shown in Fig. 1. As to sub-step 1), to provide a 
more efficient algorithm, from an initial number of identified parameters, the 

10 most significant parameters, i.e. those significantly affecting each property, 
are identified to generate a reduced number n of significant parameters, 
especially where the initial number of identified parameters is greater than 
8, as mentioned before. For that purpose, a statistical analysis algorithm 
can be used, which is based on parameter correlation calculations using 

15 parameter and property experimental data provided at prior step 42. Having 
obtained data related to / experimental runs involving an initial number p 

of parameters and a number k properties Y n each correlation factor 

contained in the correlation matrix is retained as significant whenever it 
complies with a predetermined condition in the following form: 
20 a < p {J < b or -c<p 0 <~d (1) 

wherein a,b,c and d are predetermined limit values, typically set as 
follows: 

0.5 < p & < 0.95 or -0.95 < p g < -0.5. (2) 

The parameters associated with the retained correlation factors form 
25 the reduced set of n parameters. 

It can be also shown that a minimum number / of runs at least equal 
to n + l is required to obtain reliable parameters estimation. Then, 
parameter interactions, that are in the form X i X j with i* j and which are 

significant, can be identified using the above relations (1), with the 
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suggested specific ranges given in (2). The values for X t from the / 

experimental runs are combined with the retained correlation factors p i} to 

form a final matrix W , with each element of the first column being equal to 
unity for the purpose of following sub-step 4). As to sub-step 3), it consists 
5 of establishing, for each property Y } , a best model in terms of retained 

parameters and parameter interactions. A standard variance analysis is 
carried out to confirm relevancy of all parameter coefficients and parameter 
interaction coefficients, and to select by successive variance analysis 
operations through the use of modules 24, 20 and 22, a suitable model 

10 amongst different predetermined models of upgraded degrees, whenever 
difference in performance between a given model of degree r and a 
following model of degree r+1 is found to be not significant. The resulting 
best model is taken along with matrix JFand property experimental data in 
matrix Y , as inputs for following sub-step 4) aimed at generating property 

15 behavior relations for each property Y r A matrix C of coefficient values is 

given by the matrix: 

C = (w T w)~ l xW T xY (8) 
having a dimension of [m,k], wherein m = « + / + !, t being the number of 
parameters interactions X i X j . Hence, estimated property values are given 
20 by: 



Ye =C T X = 



iW*(*.>. 



(9) 



.tX^^XfX j) 

A following step 46 as shown in Fig. 2 aimed at generating a goal 
function is carried out by the module 28 shown in Fig. 1, from the set of k 
property weights Wj produced at step 40, from the property behavior 
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relations produced at step 44 and from the specified goal values for the 
properties Y } . The basic goal vector can be expressed as follows: 

= (10) 
wherein O x is the specified goal values for the properties Y { , with / = . 
5 The goal function to be minimized may be expressed as follows: 

G(Xl,...Xn) = g T *g = t^( Y i-°i) 

which goal function is expressed in terms of property weighted deviations 
between estimated values Ye } for the properties Y } and corresponding 

specified goal values O t for the same properties Y } . A next step 48 as 

10 shown in Fig. 2 therefore consists of minimizing the goal function G to 
obtain a set of optimal parameter values for the parameters X { , which step 

48 is performed by module 30 shown in Fig. 1. Optimization step 48 
generally can consider constraints on the parameter values in the form of 
one or more ranges, typically in a form {a n b i ) wherein a^X^ b n within 

15 which optimal parameter values shall be found, according to the type of 
each variable (i.e. discrete such as binary values, or continuous). 

The "G" goal function is determinated by experimentation. The 
optimization of the "G" function is a step by step procedure. The first step 
is to obtain the behavior laws with the best fit between the experimental 
20 data and their corresponding ideal value factor. 

The second step, the optimization is based on a initial point. 

x o ={x:,xi 9 x° n \ 

and 

X M =X*-aHg\ xmxk 

25 where 

g~ v G = gradientG 
H= the Hessian of G 
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And based on the following goal function, we use the dimension reduction 
method by successive iterations 

1=1 

5 

these iterations passed by 

if / = 1 Jfc — 1 

10 and 

Now the goal function can be: 

15 G(X ] ,...,XJ=G{X„..^ ) 

We observe a perfect overlap between the two goal functions and on 
the stationary point the goal function will be; 



20 



3X; 



aofr), dG(x)df k (x^x n J 

dX ( 8X n dX t 



These equations supply the maxima and minima of the goal 

f lc (X ] ,...,X /t _ l ) 

functions including the maxima and minima from the starting goal function. 
This mathematical approach induces a reduction of the dimension of 
25 4he variables, consequently we pass from "n" variables to "n-1" variables. 



In the actual case, we start with the most important variables tror 
behavior laws with the highest weight values of the factor. 
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This approach is known under the name of network optimization, in 
this case the network nodes are built by the optimal values of the variable 
by Redressing order of the factor's rank. 

After wit^rative optimization step 48 is completed, although the set 
5 of optimal parameterValiJes ^0,. obtained can generally be considered as 
the solution to recommend/that solution is preferably evaluated amongst 
other alternative solutions by following steps 50 and 52 as shown in Fig. 2. 
At step 50, an experimental run of the process is carried out using the 
obtained set of optimal parameter values\to obtain experimental values for 
10 the k properties Y r The optimal property values X0 { and associated 

experimental property values are then evaluated at step 52 to obtain 
ranking thereof amongst a number m of other alteh^ative solutions, which 
may have been selected from knowledge base 12 shown in Fig. 1. This 
evaluation is preferably performed by a complete AHP prhcess algorithm, 
15 using the set of k property weights Wj as previously obtained through 

step 40. 

CONCRETE APPLICATION 
An example illustrating an application of the method according to the 
present invention in the pharmaceutical field will now be described. 

20 Formulation and production process for enalapril maleate tablets 

were optimized in order to provide a drug product with satisfactory 
biological performance as well as stability when packaged and stored 
under ICH (International Conference on Harmonization) conditions. Three 
(3) independent formulation and process parameters (n = 3) were identified 

25 as having an impact on the stability of the drug product: 1) the degree of 
drug neutralization during granulation {X x ) ; 2) the manufacturing 
technology ( X 2 ) ; and 3) The drug-to-excipient ratio in the formulation, i.e. 
dose strength (X 2 ) . 

As to the degree of drug neutralization during granulation (X,), it 

30 was classified as either complete, partial or no neutralization. In the case 
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of complete neutralization, the drug and the alkaline agent were both added 
to the granulation fluid, i.e. water. Therefore, the alkaline agent neutralized 
the drug prior to its addition to the powder blend for the granulation 
procedure. In partial neutralization, both the drug and the alkaline agent 
5 were added to the powder blend, blended and water added as the 
granulation fluid for the granulation procedure. When water and/or alkaline 
agent were not added to the formulation, the drug was not neutralized. The 
level of water added as well as the drug-to-alkaline agent ratio were kept 
constant for all of the formulations. The level of the alkaline agent was 
1 0 determined by the stoichiometry of the reaction. 

The manufacturing technology (X 2 ) was either wet granulation 

(X 2 =0) or direct compression (X 2 =1 ). These two technologies are used 
worldwide for the manufacturing of probably more than 90% of all of the 
solid oral dosage forms. In the wet granulation technology, the drug and 

15 other functional materials added to impart good processing attributes to the 
drug, often called excipients, are first blended together and agglomerated 
into larger particles by the addition of a granulating fluid. The role of the 
granulating fluid is to promote the development of adhesive forces between 
the materials required for the agglomeration process. After granulation, the 

20 granulating fluid is removed by drying. When a direct compression 
approach is selected as a manufacturing method, the drug is first blended 
with the excipients and tablets produced without the use of a granulating 
fluid. 

As to the dose strength (^ 3 ), four doses of the product were 
25 developed, which were obtained by using two formulations with different 
drug-to-excipient ratios (continuous parameter values) compressed at 
different tablet weights. 

A total of nine (9) experimental runs involving different formulations 
based on a combination of the three parameters were prepared, as shown 
30 in Table 1 . 
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Run 


*x 


* 2 


*3 


1 


1.73 


1 


3.33 


2 


1.73 


0 


3.33 


3 


1.73 


0 


3.33 


4 


0 


0 


3.33 


5 


0 


1 


3.33 


6 


6.9 


1 


13.3 


7 


6.9 


0 


13.3 


8 


6.9 


0 


13.3 


9 


0 


0 


13.3 



TABLE 1 



[he nine formulations covered all of the six (6) possible 
5 combinations for the wet granulation technology and three (3) combinations 
of direction compression. Tablets were manufactured by using enalapril 
maleate with USP/NF and EP excipients. In the direct compression 
technology, there is^iqt a sufficient amount of moisture to dissolve all the 
drug and alkaline agenkand provide for any significant neutralization 
10 reaction. However, excipienfcs do contain a certain level of adsorbed free 
moisture capable of creating a microenvironment where small quantities of 
, the drug and alkaline agent can be^ssolved and become available for the 

^) neutralization reaction. These phenorhena could be responsible of the 
appearance of physical as well as chemical stability problems and where 
15 taken into account by evaluating three (3) Emulation combinations. The 
nine (9) formulation combinations where prepataj and the tablets were 
stored in opened containers at 25°C/60%RH and\o°C/75%RH for a 2- 
week period. These open container studies are typicallXconducted during 
the early formulation development phases of a product to purposely 
20 accelerate physical and chemical changes in formulations in or^er to select 
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the lead candidate, i.e., the formulation with the best stability profile. After 
the 2-week time period, the tablets were removed from the environmental 
chambers and sent to the analytical department for their performance 
evaluationXThe performance of the formulations was determined by 
5 measuring teiv(^=10) properties as a function of time and temperature, 
which properties wfere selected as follows, according to a hierarchical tree 
comprising properties and sub-properties: 

- Y1 1, Y12, Y13: % drug dissolved at 5, 15, and 30 min. (sub-properties of 
10 Y1); \ 

/ - Y2: % of cyclization product at time zero; 

/ - Y31, Y32: % cyclization product affisr 2 weeks at 25°C/60%RH and at 
40°C/75%RH (sub-properties of Y3); \ 

- Y4: differential between theoretical and actual assay in mg at time zero; 
15 - Y51 p Y52: differential between theoretical a ruL actual assay in mg after 2 

weeks at 25°C/60%RH and at 40°C/75%RH (suXoroperties of Y5); 

- Y6: % hydrolytic product after 2 weeks at 40°C/75%k 

Applying the AHP process with the standard scale for these 
20 properties, the decision matrixes given in Table 2 for the properties and in 
Tables 3, 4 and 5 for the sub-properties were built. 



Goal 


Y1 


Y2 


Y3 


Y4 


Y5 


Y6 


Y1 


1.0 


(5.0) 


(5.0) 


3.0 


1.0 


(5.0) 


Y2 




1.0 


(5.0) 


1.0 


(3.0) 


(7.0) 


Y3 






1.0 


7.0 


1.0 


1.0 


Y4 








1.0 


(7.0) 


(7.0) 


Y5 










1.0 


(1.0) 


Y6 












1.0 



TABLE 2 



WO 01/18668 




PCT/CAOO/00998 



Y1 


Y11 


Y12 


Y13 


Y11 


1.0 


7.0 


9.0 


Y12 




1.0 


3.0 


Y13 






1.0 



TABLE 3 



U 
U3 



m 



CP 
u 



s 

□ 



Y3 


Y31 


Y32 


Y31 


1.0 


(5.0) 


Y32 




1.0 


TABLE 4 


Y5 


Y51 


Y52 


Y51 


1.0 


(5.0) 


Y52 




1.0 



10 



TABLE 5 

From the decision matrixes, the following weight values for the A: =10 
properties/sub-properties are given in Table 6, the sum of the weights 
being equal to unity. 



(Y11) 


w 2 
(Y12) 


w 3 
(Y13) 


W 4 

(Y2) 


w 5 
(Y31) 


W 6 

(Y32) 


W 7 

(Y4) 


(Y51) 


w 9 
(Y52) 


W 10 

(Y6) 


0.060 


0.011 


0.005 


0.098 


0.046 


0.230 


0.039 


0.034 


0.171 


0.305 



TABLE 6 
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Experimental property data that were obtained from nine (9) runs of 
the process using the selected nine (9) combinations of parameter values 
of Table 1, are given in Table 7. 



□ i in 

r\un 


V1 1 

y 1 1 


V1 0 


T IO 


T Z 


T O I 


T 


T *T 


T \J I 




Yfi 


1 


90.40 


102.7 


102.1 


0.88 


1.210 


15.06 


0.080 


(0.01) 


0.660 


4.260 


2 


91.20 


95.10 


95.80 


0.63 


1.070 


3.510 


0.120 


0.020 


0.340 


5.350 


3 


86.80 


100.4 


101.6 


0.90 


0.850 


0.720 


0.080 


0.030 


0.110 


3.670 


4 


94.00 


96.70 


97.20 


0.85 


1.340 


16.88 


0.300 


0.160 


0.780 


2.590 


5 


75.40 


106.2 


107.1 


0.83 


1.350 


22.54 


(0.08) 


0.040 


0.960 


3.570 


6 


94.70 


98.20 


98.30 


0.54 


1.110 


8.600 


0.380 


0.270 


2.200 


6.730 


7 


90.20 


98.70 


98.40 


0.64 


1.670 


5.980 


0.310 


0.290 


1.750 


7.710 


8 


51.70 


100.4 


100.8 


0.50 


0.75 


0.750 


(0.04) 


(0.14) 


0.310 


3.580 


9 


84.90 


95.30 


95.20 


0.55 


0.11 


0.110 


0.930 


0.740 


2.210 


1.450 



5 

TABLE 7 

Since n = 3 < 8, the parameter reduction step is not required for the 
purpose of the instant case. As to the statistical analysis of parameters 
10 interaction, since a correlation factor /? 13 =0.7013 for the X X X 3 interaction 
was calculated, that interaction can be considered as significant since the 
condition 0.5 < p l3 < 0.95 is satisfied. The following property behavior 
relations were established: 

15 Y11 =81.916+ 4.56 X x + 4.074 X 2 + 0.224 Z 3 - 0A23X,X 3 ; 

Y12= 101 .93-1. 45 JT, + 3.81 X 2 - 0.51 X 3 +0.14*,X,; 

Y13= 102.16 -14.8 X t +3.46Z 2 -Q.52X 3 +0A4X,X 3 ; 

Y2 =0.92-0.025*, +0.025^ -0.028* 3 +0.0018*,*,; 

Y31 = 1.42-0.23*, + 0.057 * 2 -0.03* 3 +0.019*,* 3 ; (18) 
20 Y32= 17.18-8.78*, +8.15* 2 -0.46* 3 +0.56*,* 3 ; 
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Y4 =-0.15 + 0.0193 *, -0.022 X 2 + 0.08* 3 - 0.0091 *,* 3 ; 
Y51 =-0.135-0.028*, + 0.0385 * 2 + 0.066* 3 - 0.0045 *,* 3 ; 
Y52 = 0.00089 - 0.256*, + 0.63 X 2 + 0.166* 3 + 0.008 *,* 3 ; 
Y6 =3.24 + 0.9*, +0.57*, - 0.13* 3 - 0.02 *,* 3 ; 

5 

The specified goal values for the properties as given in Table 8 were 
used to establish the goal function that was minimized to generate the 
following set of optimal parameters: 

10 *, =3.39 

* 2 = 0 (wet granulation) (1 9) 

* 3 = 7.46 

The associated experimental property values are given in Table 9. 



15 



(Y11) 


o 2 

(Y12) 


(Y13) 


o A 

(Y2) 


o 5 

(Y31) 


o 6 

(Y32) 


(Y4) 


0 8 

(Y51) 


o 9 

(Y52) 


0,0 

(Y6) 


100.00 


100.00 


85.00 


0.500 


0.500 


0.500 


0.050 


0.050 


0.050 


1.00 


TABLE 8 


Y11 


Y12 


Y13 


Y2 


Y31 


Y32 


Y4 


Y51 


Y52 


Y6 


97 


104 


106 


0.49 


0.58 


0.73 


0.031 


0.008 


0.047 


<0.05% 



TABLE 9 

20 

Applying the method for the particular case where only the minimum 
four + 1 = 3 + 1 = 4) experimental runs required were used, runs 1, 3, 6 
and 9 were selected to provide the parameter and property data as given in 
Table 7. As to the statistical analysis of parameters interaction, since a 
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correlation factor p n =0.332 for the *,* 3 interaction was calculated, that 

interaction cannot be considered as significant since 0.5 < p l3 < 0.95 is not 

satisfied. The following property behavior relations were established: 
Y11 =85.36 + 0.99 *, + 3.60* 2 -0.34* 3 ; 

5 Y12= 101.90 + 0.086*, + 2.30* 2 -0.50* 3 ; 

Y13 = 102.87 + 0.38 *, + 0.506 * 2 - 0.58 * 3 ; 

Y2 = 1 .01 3 + 0.001 *, - 0.01 9 X 2 - 0.034 * 3 ; 

Y31 = 0.88 - 0.038 *, + 0.36 X 2 + 0.0095 A", ; (20) 
Y32 = 2.88 - 2.425*, + 14.33*, + 0.61 X,; 
10 Y4 = -0.02 - 0.080*, - 0.0007 X 2 + 0.071 X 3 ; 
Y51 = -0. 0568 - 0.065*, - 0.02 X 2 + 0.06 X 3 ; 
Y52 = -0.4 - 0.08*, + 0.55 * 2 + 0.2 * 3 ; 
Y6 =2.84 + 0.68*, +0.59* 2 - 0.1 * 3 . 

15 The same specified goal values for the properties as given in Table 

8 were used to establish the goal function that was minimized to generate 
the following set of optimal parameters: 
*, = 3.32 

* 2 = 0 (wet granulation) (21 ) 

20 X, = 7.09 

The associated experimental property values are given in Table 10. 



Y11 


Y12 


Y13 


Y2 


Y31 


Y32 


Y4 


Y51 


Y52 


Y6 


78 


103 


106 


0.51 


0.55 


0.73 


0.034 


0.04 


0.029 


<0.05% 



TABLE 10 1 

25 
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Comparing the set of parameter values given at (20) with the former 
set obtained from all nine (9) experimental runs given at (19), it can be 
noted that both sets are very similar. Actually, from a pharmaceutical 
standpoint, they could almost be considered as identical. 



